Coherence Protection by the Quantum Zeno Effect and Non-Holonomic Control 

in a Rydberg Rubidium isotope 
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The protection of the coherence of open quantum systems against the influence of their envi- 
ronment is a very topical issue. A scheme is proposed here which protects a general quantum 
system from the action of a set of arbitrary uncontrolled unitary evolutions. This method draws 
its inspiration from ideas of standard error-correction (ancilla adding, coding and decoding) and 
the Quantum Zeno Effect. A demonstration of our method on a simple atomic system, namely a 
Rubidium isotope, is proposed. 

PACS numbers: 03.67.Pp Quantum error correction and other methods for protection against decoherence. 
03.65.Fd Algebraic methods. 32.80.Qk Coherent control of atomic interactions with photons. 



I. INTRODUCTION 

The uncontrollable interaction of an open quantum 
system with its environment leads to complete loss of the 
information initially stored in its quantum state. This 
phenomenon is commonly referred to as "loss of coher- 
ence". The question of how it is possible to avoid the 
negative influence of this process is one of the most in- 
teresting issues in modern quantum mechanics, and con- 
cerns many different fields of physics, in particular the 
domains of quantum information and computation. 

In the context of quantum information, the effects of 
interactions with the environment, known as "quantum 
errors" , may render information storage and processing 
unreliable 0,0- Since Shor's demonstration that error- 
correcting schemes exist in quantum computation 0, a 
general framework of error-correction has been built upon 
the formalism of quantum operations. The main con- 
tributions concern quantum codes 0], and particularly 
the class of stabilizer codes 0, 0; other strategies de- 
veloped suggest the use of " noiseless quantum codes" or 
" decoherence- free subspaces" 0,0, 13- All these methods 
usually demand that errors act independently on differ- 
ent qubits (the independent error model), and make use 
of the symmetry properties associated with these require- 
ments. This implies that the set of errors to be corrected 
hence is restricted to a special subgroup, called the Clif- 
ford group. In this paper, we present a protection method 
which draws its inspiration from the ideas of standard 
error-correction and the Quantum Zeno Effect, and re- 
quires no specific symmetry of the errors. Moreover, we 
suggest its physical implementation in an arbitrary quan- 
tum system, and show how it works for the example of a 
Rubidium isotope. 

The phenomenon known as the Quantum Zeno Effect 
(QZE) takes place in a system which is subject to fre- 
quent measurements projecting it onto its (necessarily 



known) initial state: if the time interval between two 
projections is small enough the evolution of the system 
is nearly "frozen". This effect, and its inverse (the anti- 
Zeno effect), have been widely investigated theoretically 
[Iclimri2Lll3| as well as experimentally Gener- 
alizations have been proposed which employ incomplete 
measurements 16] : in this setting, the Hilbert space is 
split into " Zeno subspaces" (degenerate multidimensional 
eigenspaces of the measured observable), and the state 
vector of the system is compelled by frequent measure- 
ments of the physical observable to remain in its initial 
Zeno subspace. The dynamics of the system in the Zeno 
subspaces has also been studied in different specific situ- 
ations 

Employing these ideas, enriched by standard tech- 
niques from coding theory [l8j . we have previously pro- 
posed an information protection scheme |l9| in Zanardi's 
spirit 20], except that we do not make any symmetry 
assumption on the unitary errors we consider. We form 
a compound system S which comprises the information 
system X to be protected and an auxiliary system A 
(called ancilla). We then apply a controlled unitary oper- 
ation C (the coding matrix) which encodes the informa- 
tion, initially stored in X, in an entangled state of X and 
A. After a short time interval, during which infinitesimal 
errors may have occurred, we apply the unitary trans- 
formation C _1 (the inverse to the previous step), which 
decodes information. Finally, we measure the ancilla to 
get rid of the infinitesimal changes that may have been 
caused by errors. Whereas in classical error-correction 
theory, the ancilla contains information about the er- 
rors allowing them to be corrected, in our QZE-based 
approach, the quantum state of the ancilla resulting from 
an elementary (coding-errors-decoding) sequence is close 
to its initial state, so that the measurement of the ancilla 
brings it back to its initial state with a probability of 
nearly 1. The key point of our method is to entangle the 
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initial state of the ancilla with the state of the system in 
such a way that the detection of the ancilla in its initial 
state implies that the system is also in its initial state. 
This is achieved through the coding procedure which has 
to ensure that after the exposition of the system to the 
action of errors, the initial state of the ancilla remains en- 
tangled only with the initial state of the system, whereas 
any other entanglement remains small during the time 
interval between consecutive measurements. 

The coding procedure involves a rather complex uni- 
tary transformation in the Hilbert space of the com- 
pound system: performing such coding operations is a 
non-trivial quantum control problem in itself. However, 
one can achieve this objective by adapting the idea of 
Non-Holonomic Control which we have previously pre- 
sented [2l|. Specific algorithms which allow us to deter- 
mine and physically implement the coding matrix C have 
been constructed. 

In this paper, we provide a comprehensive presentation 
of our theoretical scheme, including algorithmic aspects 
which were not dealt with in 

EH- 

Moreover, we propose 
a " demonstrative" application of our technique to a phys- 
ical system : more precisely, we show how our method 
can protect one qubit of information stored in the spin 
variable of the quantum state of a single Rubidium atom, 
the orbital variable playing the role of the ancilla, against 
a given set of error inducing Hamiltonians. Here we de- 
scribe a realistic experimental setting which achieves the 
different steps of our scheme through the application of a 
sequence of laser pulses and culminates in a measurement 
involving spontaneous emission. When dealing with en- 
sembles of atoms, as usually done in current cold atom 
experiments, experimental drawbacks arise due to dipo- 
lar interactions which forbid the actual implementation 
of our application. Yet, even though not completely satis- 
factory from an experimental point of view, the example 
we propose here shows the general scope of our method 
as well as its physical operationally. 

The paper is organized as follows. In Sec. II, we present 
a multidimensional generalization of the QZE and its 
application to the protection of information contained 
in compound systems. In Sec. Ill, we present the algo- 
rithms which enable us to calculate the code space and 
physically implement the coding matrix through the non 
holonomic control technique. In Sec. IV, we focus on the 
application of our method to a Rubidium isotope. In Ap- 
pendix A, we present the explicit derivation of the code 
subspace. 



II. MULTIDIMENSIONAL ZENO EFFECT AND 
COHERENCE PROTECTION 



We start this section by the geometric presentation of 
a multidimensional QZE which allows us to protect an 
arbitrary subspace of the Hilbert space against the ac- 
tion of a set of given interaction Hamiltonians. In the 
second part of this section, we take advantage of this 



phenomenon to protect an information-carrying subsys- 
tem of a compound quantum system from the influence 
of some uncontrolled error-inducing external fields. 

Consider a quantum system <S, whose ^-dimensional 
Hilbert space is denoted by TL and whose time-dependent 
Hamiltonian has the form 



M 



H{r) = E fm(r)E„ 



(1) 



where 



l ) m=l, 



are M given independent Hermi- 



tian matrices on H. and {fm(T)} m=1 M are M unknown 

functions of time. The Hamiltonian H(t) accounts for 
the errors we want to get rid of. Note that the unper- 
turbed part of the Hamiltonian is assumed to be zero 
(or proportional to the identity so that one can set it to 
zero). The standard QZE E3 El El El implies that we 
can nearly "freeze" the evolution of the system by mea- 
suring it frequently enough in its (known) initial state ; 
in other words, this effect allows us to protect the one- 
dimensional subspace spanned by the initial state of the 
system from the influence of the error-inducing Hamil- 
tonian Q. In what follows, we generalize this effect so 
as to protect an arbitrary multidimensional subspace C 
from H(t). 

Any vector \tp) of C evolves according to the operator 



U(t,t Q )=T jexp 



H{t)cLt 



where T denotes time-ordering, and where we set h = 1. 
For the QZE to hold, we shall only consider evolution 
in short time periods, whose duration T is so short that 
the corresponding action of the M components of the 



Hamiltonian is small, 
We can thus expand 



ft +T f™( T ) dT 



< 1. 



M 



U(t + T,t) = U, 



in} 



i-t+T 



fm{T)d,T E„ 



This implies that after a Zeno interval T, the initial state 
\ip) is transformed into \ip e ) = \ip) + \Stp e ) where \Stp e ) ~ 
-'Em=i £ ™-^™ W witn £ ™ = (/ /m(r)dr) . Note that, 
strictly speaking, the operator Ui n f introduced in Eq.@ 
is not unitary : nevertheless, the non-unitary part, due 
to the truncation of the time-development of the evo- 
lution operator, is of second order in time and is thus 
negligible in the Zeno limit (T — > 0). Moreover, as we ex- 
clusively consider finite dimensional systems interacting 
with classical external fields, the approximation Eq.(2) 
holds, without raising any mathematical problem. But it 
is worth emphasizing that this is no longer the case when 
dealing with systems of infinitely large Hilbert space (for 
example, see [2^1. 
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Let us assume that we are physically able to perform 
a measurement-induced projection onto C in the system 
S (see below the discussion of such projections for com- 
pound systems comprising an information subsystem and 
an ancilla). Now if we just follow the standard QZE pro- 
cedure and merely project the state vector \tp e ) , resulting 
from the infinitesimal evolution of the initial state \ip) , 
onto C, we obtain a vector \ip p }, which (a priori) differs 
from l^) (see Fig^,). This is due to the fact that usually 

the operators < E m > do not act orthogonally on 

I J m=l,...,M 

C, which means that the vectors E m and thus the in- 
crement vector \Sipe) itself are not orthogonal to C. It is 
a well-known manifestation of the standard Zeno Effect : 
the system is compelled to remain in a " Zeno subspace" 
(which corresponds here to C) in which it presents a re- 
maining dynamics called "Zeno dynamics" |17| ; in our 
case, this dynamics is unknown and thus threatens the 
information stored in C. Therefore, we see that the stan- 
dard Zeno strategy does not suffice to protect a multi- 
dimensional subspace : we have to adapt it using some 
ideas of coding theory. 

To this end, we assume a unitary matrix C acting on TL, 
which we call the coding matrix, such that the Hermitian 



operators 



{Em) 

I J m=l,. 



,M 



act orthogonally on the sub- 



space C = CC, which we call the code space. Let us de- 
note by / > 1 the dimension of C and by {|7i)} i=1 / one 

of its orthonormal bases ; ||7i) = C |7i)| will de- 

note one of the orthonormal bases of C, the state vectors 
|7i) being called the codewords. For any pair (1%) , |7 t )) 

of codewords and any operator E m £ < E m > 

I J m— 1 

have, by the definitions of C and C 



we 



,M 



(7t|7s) = fist (orthonormality condition), (3) 
7t E m 7s) = (orthogonality of the errors). (4) 



Equivalently, for any pair Qij)) , \x)) of vectors of C and 
for any operator E m £ < E m > 

I J m=l, 



M 



&E m C 1>)=0 



(5) 



In particular, for any pair (|7s) , |7t)) of basis vectors of 
C and for any operator E m £ < E m > 

I J m=l, 



li 



&E m C 7s ) =0 



If we apply the coding matrix C to the initial state vec- 
tor before exposing it to the action of the Hamil- 

tonian QJ, we obtain the new vector if}\ — C\ip) £ C 

(FigEJ)l,2) which is transformed after a Zeno interval T 



into 



U in f i> 



where 



E M c J? 
m=l £ m&n 



Em=i £ mE m C\A (Fig|n>3) 



Decoding 



ip e j yields the vector \ip' e ) = C 



■\5tf e ) where \5^ e ) ~ -i^ =1 e m &E m C \if>). From 
Eq.© it can be seen that for any vector \\) £ C, 

(xlfi^'e) = ~i E m =i £ m (xl & E m C = which means 
that \Sip' e ) is orthogonal to C (Fig|TJ)4). A measurement- 
induced projection onto C finally recovers the initial vec- 
tor with a probability very close to 1 (the error prob- 
ability is proportional to T 2 ). If the (coding-decoding- 
projection) sequence is frequently repeated, any vector 
\ip) of the subspace C can thus be protected from the 
Hamiltonian for as long as needed. We stress that the 
role of projective measurements consists both in confining 
the system in C (as in the standard Quantum Zeno Ef- 
fect) and in clearing out the erroneous component which 
has been made orthogonal to C through coding and de- 
coding. 

Let us note that a more general version of condi- 
tions Q can be considered. Indeed, if for any pair of 
codewords (|7 S ) , \jt)) of C and any error Hamiltonian 
E m £ j-E m j , 



7* E; 



fits£,m ) 



where St s is the Kronecker symbol and £ m a real number 
depending only on the number m of the error Hamilto- 
nian E m , the projection onto C — Span , i = 1, 1} 
of the state vector \ip' e ) = \ip) - !^ =1 £ ra £ ra \ip), ob- 
tained after a (coding-decoding) sequence, yields 



He W e ) = 



M 
m—1 



where H 



c 



Et=i \lt) (T*l ; if we denote by \i>) = 
a s 1 7s) the decomposition of the initial information 
state vector, HcE m \ip) has the form 



u c E m |v) = J2 a ° It*) (a 



C^E m C 



s,t=l 

I 

= £m X! a ^ ts IT* 
s,t=l 

= £m \lp) , 



(6) which finally leads to lie W e ) — fl 



* Em=l £ m£,n 



In other words, the errors E m just introduce a global 
phase factor in front of the initial information state vec- 
tor, but leaves its coherence intact. Obviously, the cor- 
rection conditions l@J are obtained as a particular case 
of the above conditions, setting £ m = for all to. Yet, 
though less general, they will be employed in the rest of 
the paper for the sake of simplicity. 
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The multidimensional generalization of the QZE we 
have just described allows us to protect any subspace 
C of a Hilbert space H against Hamiltonians of the form 
JTJ provided the projection onto C is physically achievable 
and the coding matrix C exists. This result is very useful 
in the context of information protection as we will show 
in the following paragraphs. 



Consider an information system X of Hilbert space Hi 
and dimensionality /. This system is subjected to a set 

of M error-inducing Hamiltonians \ E m \ which, 

I ) m=l,...,M 

for instance, represent interactions of the system with M 
uncontrolled external fields f m (t): we want to get rid 
of this external influence which is likely to result in the 
loss of the information stored in the initial state vec- 
tor = Yli=i c i\ v i)t where {K)}j=i,...,i denotes an 
orthonormal basis of Hi. To this end, we will use the 
multidimensional Zeno Effect. As the multidimensional 
QZE can only protect a subspace of the whole Hilbert 
space, we first have to add an A-dimensional auxiliary 
system A (called ancilla) to our system X, so that the 
information is transferred from Hi into an /-dimensional 
subspace C of the (N = I x v4)-dimensional Hilbert space 
H = Hi®Ha of the compound system S = 1 ® A (let us 
note that this ancilla adding procedure is quite standard 
in quantum error-correction [l| and corresponds to the 
redundancy used in classical coding methods). Further- 
more, we shall suppose that all the state vectors of the 
different Hilbert spaces Hi, Ha and hence H are degen- 
erate in energy so that the unperturbed part H of the 
Hamiltonian can be set to zero as in the first part of 
this section: the subspace C and the information it car- 
ries can thus be protected through the multidimensional 
QZE (in Sec. IV, on the example of Rb, we shall see that 
the multidimensional QZE may also be used even though 
Hq is not zero, provided Hq and the errors have some 
convenient properties). Note that A and X need not be 
"physically separate" systems, but only have to possess 
independent Hilbert spaces Ha and Hi- For example, 
in Sec. IV, we shall consider the Rubidium atom as the 
compound of two independent subsystems, namely its 
spin (which plays the role of X) and orbital part (which 
plays the role of .4.). Doing so, we shall use the terms 
'factorized' and 'entangled' in a generalized manner to 
designate states obtained as a direct product of the spin 
and the orbital parts, and linear combinations of such 
states, respectively. 

Let us now return 
sider the simple case 
in the pure state \a). 
ried by \ipi) £ Hi is 
torized state 1^) = \ipi 



to our problem and first con- 
in which the ancilla is initially 
The information previously car- 
then transferred into the fac- 

> ® la) = ELi c » \ v i) ® M = 



Si=i c i which belongs to the tensor product subspace 



{|7<> = Wi) <8> |«)} i=w 



C = Hi®Span\\a)] = Span 

In other words, the initial density matrix of the com- 
pound system S is p = {ipi\) <8> (|a) After 



coding (through the matrix C) it reads p = C^pC ; 
at the end of the action of the errors it is transformed 
into p e = U- n fC^pCUi n f ; finally it takes the form 

p e = CU^.C^pCUinfC^ after decoding. In this setting, 
the projection onto C can be simply achieved by measur- 
ing the ancilla in its initial state \a). As T is very short, 
the state of the ancilla evolves just a little within a Zeno 
interval, such that the probability of detecting it in its 
initial state \a), and thus of projecting the state of the 
compound system onto C is very close to 1. After projec- 
tion, we trace out the ancilla to obtain the final reduced 
density matrix p? I = (a\ CUj n ^C^pCUi n fC^ \a) for the 
information system X; in the same way, one can calcu- 
late the initial reduced density matrix is pi = \ipi) {ipi \ . 
The variation 5pi = pfj — pj of the information-space den- 
sity matrix during the whole process can be expressed as 
the commutator 



Spi = -i 



M . 

J fm{r)dT {a\C^E m C\a),pi 



from which we infer that pi satisfies the equation 



. dpi 

1 dt 



[he, Pi] , 



A I 



he = ^2f m {a\d^E m C\a) 



m— 1 

where h e is an effective Hamiltonian which is determined 
by the error-inducing Hamiltonians transformed by the 
coding and decoding and projected onto the initial state 
of the ancilla. From Eq.© one can infer that h e = 
and hence pi remains constant in time: as long as we re- 
peat the coding-decoding-ancilla resetting sequence, the 
information initially stored in X is protected. 

It is not always feasible to directly measure the an- 
cilla independently from the information system ; in 
other words, it is sometimes impossible to perform a 
projection onto disentangled subspaces of H of the form 
Hi® Span \\a)] : in some cases, as for the Rb atom 
(Sec. IV), one can only project onto entangled subspaces 
of the total Hilbert space H. In such a case the informa- 
tion initially stored in the vector \ipi) = X^<=i °i \ v i) 
Hi is transferred into an entangled state of X and A 
of the form |V>) = J2i=i c « |t<) where the I vectors 
(i = 1, which form an orthonormal basis of the 

information-carrying subspace C, are not factorized as 
earlier but are in general entangled states. Nevertheless 
the same method as before can be used in that case to 
protect information, albeit in a different subspace C. 

To conclude this section, let us make a few re- 
marks about our method. We first emphasize that our 
technique, though inspired by quantum error-correcting 
codes U, is very different from them : indeed, in those 
schemes, the information is encoded in such a way that 
it can be corrected from the action of a set of errors 
through a syndrome measurement, followed by a (con- 
ditioned) recovery operation, depending on the result of 
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d ) Without coding -decoding 




b ) With coding-decoding 




FIG. 1: Multidimensional QZE: a) a simple projection fails to recover the initial vector, b) the sequence coding-decoding- 
projection protects the initial vector. 



the measurement; on the other hand, in our technique, in- 
formation is continuously protected by the frequent rep- 
etition of a three step cycle (coding-decoding-projective 
measurement), in which the projective measurement does 
not give any indication about which error occurred, but 
simply clears out the erroneous component of the state 
vector, which has been made orthogonal to the initial 
information-carrying subspace through coding and de- 
coding. Let us now return to conditions and J3J im- 
posed on the codewords ,i = 1, /} and make two 
points about them: 



A. We can establish a useful relation between the di- 
mension of the ancilla and the number of correctable er- 
ror Hamiltonians. The set of the I codewords can be seen 
as a collection of 21 x N = 2I 2 A real numbers on which 
2T 2 + 2MI 2 = 2I 2 (l + M) constraints, directly derived 
from Eas. (|HI4J) . are imposed. The number of free parame- 
ters must be larger than the number of constraints, hence 
we necessarily have 2I 2 A > 2I 2 (\ + M), which satisfies 

A — 1 > M. (7) 

This condition gives an upper-bound on the number of in- 
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dependent error-inducing Hamiltonians that our method 
can correct simultaneously and is called the "Hamming 
bound" . 

B. We may compare our correctability conditions Q 
with the more general conditions (see [l| p. 436) of stan- 
dard quantum error-correction 



V(|7 s ),|7 f »GC 2 , V(E fe ,E ; ) e{% ({E m })}, 

7s) = an (jtlls) (8) 



7t 



which ensure the existence of a code space that is com- 
pletely protected against the error-inducing Hamiltoni- 
ans E m . Here otui are complex numbers, and the set 
| -Em | of Hermitian operators E m generates a group 

Q d^" 1 }) °f au possible error- induced evolutions J2J). 
By | Ej ^ | E m | ^ | we denote a complete basis set of op- 
erators which spans the space of evolution operators U 
and allows one to represent any U as a linear combina- 
tion of the basis operators Ej . The variety of all linear 
combinations of E 3 includes not only all E m but also 
many other operators given by commutators of all or- 
ders in E m entering the expansion of U for long times. 
The condition © is therefore much more restrictive than 
Eq.Q. Moreover, even for just two generic matrices E m , 
the basis {^j j spans the entire Hilbert space TL, yielding 

C = 0. Only if the set |-B m | belongs to an extraspecial 
algebra restricting the error evolution operators U to a 



subgroup Q ({-^m}) Su (W) of the full unitary group 

in 7Y, a non-trivial code space C may exist. The Zeno 
effect is the only way to suppress loss of coherence if it is 
not the case. 



III. THE CODE SPACE AND THE CODING 
MATRIX 

It is sometimes possible to build the code space C ex- 
plicitly from physical considerations: Appendix A gives 
an example of a situation in which the code basis can be 
found directly. In general, however, we need an algorithm 
to calculate the code basis {|7i)} i=1 7 or, equivalently, 

the coding matrix C. We start this section by describing 
this algorithm. Then, in a second part, we show that the 
non-holonomic control technique |2l| can be employed to 
implement the coding matrix physically. We also provide 
an algorithm which achieves the appropriate control. 

Let us first make a remark which will be useful in what 
follows. Consider a vector |C) of some Hilbert space and 
a matrix E on this space. From the vector |C) we want 



to calculate a vector 



C 



C) such that (C 



C) = 0, then |C) = 



/c(A) 



C 



C ) and the function 



AE 



0. If 



depending on the c- number A, is minimal for A = 0: in- 
deed 



AE 



C|C) + A(C 



A* (C 



C) + |A| 2 (C E+E C) = 1 + |A| 2 (C 



EtE 



and as ( C 



EtE 



is A = 0. But if ( C 



C) > 0, fe(A) is minimal for |A| = 0, that 



7^ 0, we can apply the following 
iterative method: we minimize /c(A) with respect to A, 
then we set |C) = |C) + -|E |C) and take -^J= as our 
new |C) ; we repeat this sequence as long as needed: |C) 



C ) , such that ( C 



C> =0. 



finally tends to 



Let us now return to our problem and show how the 
previous remark can help us. What we want is to find 
/ vectors which meet the conditions J3J| and Q ; 
equivalently, we can say that we look for an orthonormal 
basis in which all the matrices Ej~ have their I x I upper 
left blocks equal to zero. To solve this problem, one can 
first be tempted to use standard techniques of linear al- 
gebra, in particular matrix diagonalization : however, it 
appears that these methods do not work, except in the 



trivial case when all the matrices E& have a common ker- 
nel, which is much more than what conditions © and Q 
require. So we propose to transform our initial problem 
in such a way that it can be dealt with by the iterative 
algorithm presented in the previous paragraph. Let us 
combine the / vectors |7,) into a (N x /) " supervector" 



C = 



Then let us build E = 



fill) 



\l7/> 



M 



different 



(TV x I) x (N x 7)-dimensional super-matrices E^ in the 
following way: we consider them as made of I 2 blocks of 
dimension N x N and we successively fill each of these 
blocks with the different Hamiltonians E m or the iden- 
tity matrix / or 0. To be more explicit, the first I ^ I „ 1 - ) 



7 



matrices are built by simply placing the N x N identity 
matrix in each of the blocks situated above the 



MI(I+l) 

2 

and above the diagonal. One can thus reformulate the 
conditions J3J as follows: for 1 < k < / ^ / 2 ' 1 - ) 



diagonal. In the last 1 ^LL1L1±> ones, the M operators E, 
are successively placed in each of the blocks on 



will be the starting point of the first step: we normal- 
ize this vector by imposing to each of its I components 
to have norm = j. If one of the components of |Co) is 
non normalizable, that is equals zero, we pick up a new 
random supervector | Co) as a starting point. 

Then, as in our preliminary remark, we minimize the 
function 



C 



0. 



Note that this form does not take the normalization con- 
dition into account, which will be imposed differently. 
Similarly, the conditions Q are translated into the fol- 



lowing form: for 



1 < k < 



2 



M 1(1+1) 



c 



0. 



Our initial multivectorial problem given by Ea . 13I4JI has 
thus been transformed into a simpler one which can be 
handled by the same kind of iterative algorithm as in 
our preliminary remark: we just have to find a (N x I)- 



dimensional supervector 

1(1-1) , M 1(1+1) 



C) such that for 1 < k < 



0. 



k=l 



|Co) + Ai 0) E fc |C ) 



with respect to the E c-numbers A? : actually, we sepa- 



rate the real and imaginary parts of A^ "* = 



1.0) 



and calculate the appropriate 's and (3^ 's by solving 
the set of 2E equations 

dF Co n 



dF Co 



0, 



which can be translated into the linear system 
X(|Co)).^ (0) =^(|C )) 



Let us now review our iterative algorithm in more de- 
tail. First we randomly pick a supervector | Co) which 



where K (|Co)) is a 2E x 2i?-dimensional real matrix de- 
fined by 



Re (Co 



Kij (|C » = I 



-Im ((C 
Imf/Co ' 



Re ( C 



EtE, 



ELSE; 



-E j — E 



Co 
Co 
Co 
Co 



for 1 < i < E and 1 < j < E 
for 1 < i < E and 1 + E < j < IE 
for 1 + E < i < 2E and 1 < j < E 
for 1 + E < i < 2E and 1 + E < j < 2E 



r 



D (| Co)) is a 2i?-dimensional real vector defined by 



^(|Co» = 



-Re 
Im ( / Co 



-i-E 



Co 



for 1 < i < E 
for E + 1 < i < 2E 



>(o) . 



and A is a 2£'-dimensional real vector containing the 



parameters oti °''s and ^"""s 



,(0), 



t (0) = 



M 

(0) 



Once the e-immlnTS A^ = 



(0) 



•4 0) 



found, we calculate |ACo) = Ylk^k -ft I Co) and |C ) = 
| Co) + \ |ACo). We normalize |C ) by requiring each of 
its / components to have the norm = j, and take the 
result of this operation as our new starting point |Ci). If 
one of the components of |C ) is non normalizable, that 
is equals zero, we pick up a new random supervector | Co) 
as a starting point. 

We repeat this sequence of operations as long as 
needed. Thus, at the mth step, we minimize the function 



's have been 



x(m-l) \(m-l) 
A l i A 2 > 



,A 



(m-l) 
E 



E 

E 

k=l 



|C m _ 1 ) + Ai m - 1) E fe |C 



771 — 1/ 
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by solving the real linear system 



*T(|Cm-i»-A 



(m-1) 



= ^(|C m _ 1 ». 



This yields the A£ 



(m— 1) 



s and 



AC m _i) from which we 
calculate |C^_ X ) = |C m _i) + \ |AC m _i). If possi- 
ble, we normalize |CJ n _ 1 ) and take the resulting vec- 
tor as the starting point |C m ) of the (m + l)th step 
; otherwise, we pick up a new vector |Cq) as a start- 



ing point. Finally |C m ) tends to 



C) such that Vfc € 



C 



C) = 0. 



This algorithm was numerically implemented and al- 
lowed us to exhibit new codes : we protected 2 qubits 
among 7 against the action of 31 errors (21 individual er- 
rors + 10 collective errors) and 4 qubits among 9 against 
the action of 27 individual errors |l9{|. 



The coding matrix C which allows us to transfer the 
information from the space C to the code space C, is a 
rather complex unitary operator on the Hilbert space of 
the compound system <S = T ® A. We have just shown 
how to calculate the codewords, which actually form the 
first / columns of C, but one can wonder how to imple- 
ment it physically. The question of the physical feasi- 
bility of the coding matrix C can be solved by the non- 
holonomic control technique. 

The non-holonomic control technique has been sug- 
gested by our team as a means of controlling the evo- 
lution of quantum systems |2l| . Basically, it consists in 
alternately applying two "well-chosen" perturbations V a 
and Vb to the system S we want to control during pulses 
with timings U. The total Hamiltonian H = Hq + V 
thus has a pulsed shape and alternately takes the two 
values H a = Hq + V a (during odd-numbered pulses) and 
Hf, = Hq + Vb (even- numbered pulses). The timings U 
play the role of free parameters one has to adjust in or- 
der to perform the desired control operation. To be more 
explicit, the perturbations V a and Vb must be chosen so 
that the commutators of all orders of H a = Hq + V a and 
Hb = Hq + Vb span the whole space of Hermitian ma- 
trices acting on the system we want to control: this is 
called the bracket generation condition (BGC). From the 
Campbell-Baker-Hausdorf formula, it follows that this is 
a necessary condition of controllability. It also proves 
to be sufficient in all the practical cases we dealt with. 
For that reason, we consider that we have "good con- 
trollability conditions" as soon as BGC is checked. The 
number nc of control timings depends on the problem 
to be solved. For instance, if we want to impose the ar- 
bitrary evolution U ar b on an iV-dimensional system, we 
need at least nc — N 2 timings ti, since iV 2 is the total 
number of free real parameters characterizing a N x N 
unitary matrix. We dealt with this problem of complete 
control in previous papers |2l| . and developed a general 



algorithm to find the appropriate timings ti which realize 

U (ti,*2, —jijv 2 ) = CX P (-iH a t N 2) exp f-iff&tjv 2 -i) 

. . .exp (-iHbtl) = U arb . 



We can directly apply this result to our coding prob- 
lem in the following way: first, we find the codewords 
7 i — 1) ■■■7 1} by the iterative algorithm we have pre- 
sented in the first part of this section, then we complete 
the set of I vectors {|t$) ,i = l, —,1} with (N — I) vec- 
tors {{'Jj} ,j = I + 1, to form an orthonormal basis 
of Ti, we build the coding matrix by taking the vectors 
{|7») 7 i = 1) •••! N} as columns of C, and finally we calcu- 
late the nc = N 2 appropriate timings {ti} such that 

U (*i,*2, ■■■,t N 2) = exp (-iH a t N 2^j exp i^-iHbt^^ 

. . . exp (^—iHbtij = C 

through the complete control algorithm presented in |2l| . 
Note that we assume Hq — (Sec. II), hence H a = V a and 
H b = V b 

Actually, this procedure provides a lot of useless work: 
indeed, most of the information contained in the coding 
matrix is irrelevant and the N 2 real parameters of C do 
not all have to be controlled exactly: the number nc of 
necessary control parameters {ti} is much less than N 2 . 
Let us examine this point in more detail. 

The coding matrix is characterized by the relations 
©. The problem of control thus reduces to finding nc 
timings ti , which we will formally gather in a time- vector 

t = I , such that the non-holonomic evolution 

\ 'nc / 

matrix 

U (T) = exp (-iHatn^j ex P (-*#&*rjo-lJ 

... exp (-iH a t^J 

meets conditions (0: for any pair (\j s ) 7 1 Tt ) ) i< s t<i °f 

basis vectors of C and any operator E m G \ E m \ 

I J m=l, 



It 



C/ t (?) E m U (?) 



7.= 



0. 



(9) 



The number nc of control parameters must exceed the 
number of independent constraints which is clearly ~ 
MI 2 , that is nc ^ MI 2 . The number of really necessary 
control parameters appears to be much smaller than N 2 . 
We have to design a new algorithm which achieves a par- 
tial and less expensive control of the evolution operator 
of the system. 

The algorithm we shall use to calculate the appropriate 
control timings ti mixes the iterative algorithm presented 
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at the beginning of this section and the non-holonomic 
control technique. If we introduce the (N x I)x (N x I)- 
dimensional block-diagonal matrix 



0(T) = 



fu(t) 
u(~t 



\ 



\ 



o o ■■■ u(T) J 

and the (N x J)-dimensional supervector 

/l7i 



IQ = 



V hi) 



composed of the coordinates of the / basis vectors of C, 
we can set the problem of control Eq. © in the following 
equivalent form: we look for a time-vector t such that 



Vfc, 



c 0t(T)Efc0(T) c) = o 



where the matrices 



i j fc=i,...,i 



(10) 



denote E different ma- 



trices of dimension (N x I)x(N x I) which have been in- 
troduced in the beginning of this section. In other words, 
we look for the time-vector t which sets to zero the 



test function G(T 



C 



ut(T)E fe 0(T) 



The idea of our algorithm is to take the super vector 
|Co) = U ( t o) |C), where t o is a random time-vector, 
as the starting point for an elementary step of the it- 
erative algorithm and look for the small time increment 
dto such that U (l?o + dto^j |C) follows the direction pro- 
vided by the result |Co) + |ACo) of the iterative algo- 
rithm. The repetition of this sequence finally yields 
t = T + dt + dt\ + ... which meets Ea. ljTU)) . 

Let us now describe the algorithm in more detail. 
First, we randomly pick a set of timings to.; in a "re- 
alistic range" , dictated by the system under considera- 
tion: in particular, control-pulse timings have to be much 
shorter than the typical lifetime of the system and be 
much longer than the typical response delay required by 
the experiment. Then we minimize the function 



Fc (A? 



(0) x (0) 



v (0) 



fe=l 



|C > + Ai 0) E fc |C 



as we did in the algorithm presented at the begin- 
ning of this section: we obtain the A^'s and |ACo) = 
^ fc AfeEfc |Co). At that point, we look for the small incre- 
ment dto of the time- vector t o such that 



Vfc, /c (T„) .dto J E fc U (T ) + Qt (To) E fc (™ 



Co + iACo % Co + ±AC 




(Co + iACo|Co + iACo) 



(11) 



It should be noticed that we do not consider the 
error super-matrices E& corresponding to orthonormal- 
ity conditions: in other words, we just take matrices 

{M fee [^ +1 ,£^ + -<^] into accmmt - Thuswe 

deal with MI (* +1 ) complex equations. This set of equa- 
tions can be reduced to the real linear system 

5 (T ) -dto = ^(|AC )) (12) 

where S (To) and W (|ACo)) are respectively an MI 2 x 
nc real matrix and a M/ 2 -dimensional real vector. We 
obtained Eq. (|12fl by splitting the set of M7 ^ +1 - ) complex 
equations (|ll|l into two sets of MJM+H Tea \ equations, 
and rejecting those which are trivial (0 = 0) or redun- 
dant. Even though this procedure is straightforward, the 
explicit expressions of the different elements of 5* and W 



involve many indices and are so unpleasant that we prefer 
not to reproduce them here. 

The linear system we have just found is a priori rect- 
angular (MI 2 x nc), but actually we have not fixed the 
number nc yet. Previously, we stated that nc > MI 2 : 
we could be tempted to set n c — MI 2 so as to obtain 
a square system, easily solvable by standard techniques 
of linear algebra. Yet we will proceed in a slightly dif- 
ferent way. We set nc > MI 2 , say nc — MI 2 + 5n 
where Sn is an integer of order 1. Then we randomly 
pick MI 2 timings ti among the nc which will be consid- 
ered as free parameters, whereas the other Sn ones will be 
regarded as frozen. In other words, we randomly choose 
a permutation <jq G S nc (symmetric group of order nc) 
and take the timings {ij = i CTo (o}—i MJ2 as free pa- 
rameters whereas the timings -ft' — t rrn i i \ } . are 

° L 1 °0\ l J i z—l + M I z ,nc 

frozen. This leads to new versions of Eas (|llll2fl : 
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Vfc, / C f ^ (T ) ■ dt' j E,U (T ) + Qt (To) E fe ( 



Co + ±AC 



Co 




rACo)- 



(Co + iACo|Co + iACo) 



(13) 



S(T )-d?o = ^(|AC )). 



(14) 



Eg. 114(1 is now clearly a square system. Solving Equa- 
tion l|14(l yields the M7 2 -dimensional increment dt'o 
which we complete with 6n zeros into a n^-dimensional 
vector ; by reordering timings we obtain the total 
time- vector increment dtg. Thus we have for i G 
[l,M/ 2 ], cZto,a- (i) 7^ (free parameters), whereas for 
i G [l + Af/ 2 ,n c ], rfto,<T (j) = (frozen timings). Then 
we set T" = To+ a dto where a is a convergence 
coefficient and calculate the test function G ( t ) = 



in t 



t ! for differ- 



J2 k (c ut(T)E fe u(T) 

cnt values of a e [0,1]. If we find an a% such that 
G ^ < G (To) , we take Ti = Tj 1 as our new time- 
vector, and keep the same free- varying timings: in other 
words, the permutation o\ governing the timings that 



play the role of control parameters in the second step of 
the algorithm remains the same, that is o\ = <tq. If we 
cannot find an appropriate ot\, this means we are situ- 
ated in a local minimum of G ; then we set t \ = to 
and pick a new set of free varying parameters by simply 
choosing a new permutation o\ ^ (Tq randomly. This 
rotation procedure among control parameters allows us 
to avoid possible local minima of the test function G we 
want to cancel. 



We repeat this sequence of operations as long as 
needed. At the rath step, we take the supervector 
|C TO _i) = U ( t m -i) |C) as the starting point of an 
elementary step of the iterative algorithm. We cal- 



culate |AC r 



= £A 



(m-l) 



and find M I 2 - 



dimensional variations vector dt' ' m -i of the MI 2 free pa- 
rameters (characterized by permutation cr m _i) such that 



Vfc, ( C 



^ (T m ) • dt'rnj EfeU (T m ) + 0* (T m ) E fc (T m ) ■ dt' 

C m -i + ^AC m _i Ek C m -i + iAC m _i) — (C TO _i 




(C m -i + j AC TO _i|C TO _i + iAC m _i) 



by solving the associated square linear system 

S (? m -l) ■ d?m-l = & (|AC m _!)) . 

We complete dt TO _i with 5n zeros and reorder the tim- 
ings so as to obtain dt m -\. Then we take t m = t m — 1+ 
a dt m -\. If there exists an a m such that F ^T m " 1 ^ < 

F (t m _i) we set t m = T m m as our new time- vector, 
and keep the same free parameters for the (to + l)th 
step: the permutation characterizing free- varying timings 
in the (ra + l)th step will be the same as in the rath step, 
that is <r m = <j m -\- Otherwise, we take t m — t m -i 
as our time-vector and randomly pick up MI 2 new free 
parameters among the nc timings, by choosing a new 
permutation <j m for the (to + l)th step. 

We have not said anything about the decoding so far. 



If the two Hamiltonians H a — V a and iff, = V& can be re- 
versed (note that we assume H = 0), i.e. the sign of V a 
and T4 can be reversed by altering of the control field pa- 
rameters, the implementation of the decoding matrix is 
quite easy: it amounts to reversing V a and Vfc and apply- 
ing the same control timing sequence backwards. To be 
more explicit, one starts by applying —Vfc during timing 
t nc , then —V a during i„ c _i,... , and finally — V a during 
t\ . On the contrary, if V a and Vfc cannot be reversed, one 
cannot apply this technique. We must use the general 
non-holonomic control technique, involving ./V 2 control 
parameters, to find timings which realize C~ l . 

The algorithm we have just described was numerically 
implemented and has already given satisfying numerical 
results on a realistic 7-qubit system subject to the action 
of 21 errors In the next section, we deal with another 
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real physical system which lends itself particularly well 
to a demonstration of our method. 

To conclude this section, let us emphasize that, to 
our knowledge, there is only a formal link between our 
method and the so-called "bang-bang" control schemes 
[23l |. Actually, in this kind of techniques, fast and strong 
pulses are applied which average the interaction Hamil- 
tonian between the system and its environment to zero. 
By contrast, our method employs pulses which are de- 
signed to code information, that is to transfer it into a 
proper subspace, in which errors act orthogonally : de- 
coding and measurement then allow us to recover initial 
information. 



IV. COHERENCE PROTECTION APPLIED TO 
THE RUBIDIUM ATOM 

The goal of this section is to apply our method to a 
real physical system. As we shall see below, the cho- 
sen system, a Rubidium isotope, due to its structure, 
lends itself particularly well to a straightforward imple- 
mentation of our technique and allows us to illustrate its 
different steps quite simply : to be more specific, follow- 
ing the scheme we presented in the previous sections, we 
show that it is possible to protect one qubit of informa- 
tion encoded on the two spin states of the ground level 
5s of the radioactive isotope 78 Rb against the action of 
M = 6 error- inducing Hamiltonians E m . For numerical 
calculations we considered 3 magnetic Hamiltonians 

7p 



oc L k + 2S k ,k = x,y, 



and 3 electric Hamiltonians of second order 



E 



kj 



i k.l — 



x, y, z 7 



k<v, 



In the following, we propose a detailed physical setting 
which achieves the desired protection operation : in par- 
ticular, we provide characteristic values of control fields 
and pulse timings. These different calculated parameters 
relate to a single isolated atom. As we shall see at the end 
of this section, when dealing with an ensemble of atoms, 
serious experimental drawbacks emerge which prevent us 
from actually implementing our application. Neverthe- 
less, the example considered shows the operationality of 
our method which is able, in a given physical situation, 
to provide a precise frame for its implementation. 

Before presenting the details of the proposed imple- 
mentation, let us motivate the choice of the Rubidium 
atom. Alkali atoms like Rb are very interesting for our 
purpose because of their hydrogen- like behavior. Such an 
atom is the compound of an information subsystem, i.e. 
the spin part of the wavefunction, and an ancilla, i.e. the 
orbital part of the quantum state. As we shall see, it is 
easy to increase the dimensionality of the ancilla by sim- 
ply pumping the atom towards a shell of higher orbital 
angular momentum L. 



60f,T*0M5i 

7956.9 cm-' 

1.26 //tn 
0.9865 eV 



J=5/2 



<I 



2,3 ia 5 cin-i = 310-"eV 



J=7/2 



J=5/2 



5d,rx266.2 ns \ 

1.2531 la 4 are' 
1.5536 10-x eV 



3 cm' 1 

0.000372 eV 



1.5 cm' 1 
1.5 cm' 1 



13004 cm' 1 

769 nm 
1.61 eV 



5p,T*26.51 ns V 

1.258 10-' cm' 
1.5601 10 7 eV 

12697.7 cm 1 

788 nm 
1.5743 eV 



J =3/2 
J =3/2 



237.6 cm 

0.0295 eV 



^18.8 c 



118.8 cm- 



5s 



J=1/2 



J=1/2 



FIG. 2: Spectrum of 78 Rb: The useful part of the spectrum 
of Rubidium is represented. 



We chose 78 Rb among all alkali systems because of 
its spectroscopic characterics (Fig|2J) [2J,|23. In partic- 
ular, 78 Rb has no hyperfine structure (its nuclear spin 
is 0) which ensures that the ground level 5s is degen- 
erate: this is necessary for the projection scheme as we 
shall see below. Moreover it has a long enough lifetime 
(t ~ 17.66 min) for the proposed experiment. 

Let us now review each step of our method in detail. 
As mentioned above, the information we want to pro- 
tect is initially encoded on the two spin states = 



\5s,j = ±,m.j = -i) and \v 2 ) = \5s,j = \,m 3 = ±) of 
the ground level 5s of the atom: these two states span 
the information space Hi = Span[\vi) , |^)] whose di- 
mension is in that case 1 — 2. The first step of our 
scheme consists in adding an ancilla A to the information 
system. The role of A is played by the orbital part of the 
wavefunction. In the ground state (L = 0), its dimension 
is A = 2L + 1 = 1 (roughly speaking, there is no ancilla) . 
If we want to protect one qubit of information against 
M = 6 error- inducing Hamiltonians, we have to increase 
the dimensionality of the ancilla up to A = M + 1 = 7 
(Eq.(J7|): this can be achieved by pumping the atom up 
to a shell nf (L = 3). We choose the highly excited 
Rydberg state 60/ so as to make the fine structure as 
weak as possible (the splitting for 60/ is approximately 
10 _5 cto _1 24]). We shall first consider the fine structure 
is negligible so that the N = I x A = 2x7 = 14 basis vec- 
tors of the total Hilbert space TC = Tti <S> Ha are almost 
perfectly degenerate ; the validity of this approximation 
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will be discussed at the end of this section. To be more 
specific, the pumping is done in such a way that 



^2 



I7i) 
I72) 



60/, j = -,mj 



60/, j 



In other words, using the terminology of the previous sec- 
tions, the information initially stored in Tij is transferred 
into 



C = Span 



71 



60 f, 3 



-,m,j 



60/, j 



The choice of the subspace C may appear arbitrary at 
this stage, but it will be justified later by the practical 
feasibility of the projection process onto C. Let us note 
that C is an "entangled" subspace whose basis vectors 
{l7«)}i=i 2 are general entangled states of the spin and 
orbital parts: this means (Sec. II) that the projection step 
will not consist in a simple measurement of the ancilla 
but will involve a more intricate process we shall describe 
in detail later. 

Practically, the pumping can be achieved as follows. 
One applies three lasers to the atom: the first laser 
is right polarized and slightly detuned from the transi- 
tion (5s <-> 5p) whereas the second and third lasers are 
left polarized and slightly detuned from the transitions 

^5p| and ^5di <-> 60/^j respectively. The de- 

tunings forbid real one-photon processes: the atom can 
only absorb three photons simultaneously and is thereby 
excited from the ground level 5s to the Rydberg level 60/. 
By using selection rules, one can construct the allowed 
paths represented on Fig|3] these paths only couple 
and I V2) to j 71) and I72), respectively. 



whose frequencies lor and uj' r are respectively slightly 
detuned from the two transitions (60/ <-> 5d,j = §) and 
(60/ <-> 5d, j = |) (detunings S and 5'). The character- 
istic values of these fields are 

E x , a = E' xa = 8.5 ltfV.m- 1 



Ey t a — 

<Py,a = 

TlUJr — 

s = 

huj' R = 

6' = 



E' ya = 5.2 ltfV.m- 1 
0.986324 eV = 7955.14 cm" 



-0.000010 eV 
0.986676 eV 
0.000010 eV 



-0.080654 cm' 
7958.14 cm -1 
0.080654 cm" 1 . 



The intensity of the laser beams are typically of the or- 
der of 2 10 8 W.cm~ 1 . The Raman Hamiltonian associated 
with these fields is denoted by Wr.a- The total perturba- 
tion is V a = Wz + Wr a- During even- numbered pulses 
("B" type pulses), we apply the same magnetic field as 
for A type pulses, which is experimentally convenient, 
and two sinusoidal electric laser fields 



The second step consists in encoding the information 
by the non-holonomic control technique: to impose the 
coding matrix on the system, we submit the atom to 
nc = 34 control pulses of timings {ti} i=1 34 , during 
which two different combinations of magnetic and Raman 
electric Hamiltonians are alternately applied (see FigQJ. 
To be more explicit, during odd- numbered pulses ("A" 
type pulses) we apply a constant magnetic field 



E* b (t) = Re \%, 
where ~Eu 



^ — iUJryt 



Ex b 







K.b 

yj> ' 







whose frequencies are the same as above and whose char- 
acteristics values are 



B x = 7 1CT 3 T 
B y = 8.2 1(T 3 T 
B z = -6.8 1CT 3 T 



which is associated with the Zeeman Hamiltonian Wz, 
and two sinusoidal electric laser fields 



#«(t) = Re [^e-^'l ,X(t) - Re [He 



^—lUj'jyt 



En 



K 




x.a 



E x ,b 

Ey.b 



EL 



-5.2 10 6 y.m" 



E' y b = 8.5 lO^.m" 1 



<Py, a 



2.3. 



The Raman JTamiltonian associated with these fields is 
denoted by Wr,b- The corresponding perturbation is 

Vb = Wz + Wr.b- Therefore, since the fine structure of 
the level 60/ is neglected, the unperturbed Hamiltonian 
Hq is taken to be and the total Hamiltonian has the 
form: Ha = V a during "A" pulses, Hb = Vb during "B" 
pulses. The 34 different timings have been calculated so 



13 




FIG. 3: Ancilla adding by Pumping. Photon polarization and involved sub-Zeeman levels are represented. The fine structure 
of the Rydberg level 60/ is not resolvable. 



that 

U(tl, ...,t 34 ) = e -iHBt noe -iH A t no -i e -iH A t x = Q 

meets conditions ©. At the end of the coding step the 
information is transferred into the code space C = CC, 
encoded on the codewords \ ]%) = C Kj) \ 

I J i=l,2 

As can be easily checked from Fig0]the total duration 
of a control period (~ 125ns) is approximately 10 3 times 
shorter than the lifetime of 60/ Rydberg state which is 
approximately 0.115ms as can be calculated from [24| . 
The different pulse timings range between 2.9ns and 
7.4ns, which are feasible. 

After a short time, the information stored in the sys- 
tem acquires a small erroneous component due to the 
action of the error Hamiltonians, which is orthogonal to 
the code space C. Then, we apply the decoding matrix 
C^ 1 to the atom as suggested at the end of Sec. III. We 
reverse B and the detunings 5 and 5' , and leave all the 
other values unchanged (this amounts to taking the op- 
posite of Hamiltonians Ha and Hb), and apply the same 
sequence of control pulses backwards: we start with an 
"A" pulse whose timing is t nc , then apply a "B" pulse 
during i„ c _i, etc. (see Fig|3J). The decoding step yields 
an erroneous state whose projection onto C is the initial 
information state. 

In the last step the erroneous state vector is pro- 
jected onto the subspace C to recover the initial infor- 
mation. Projection is a non-unitary process which can- 



not be achieved through a Hamiltonian process, but re- 
quires the introduction of irreversibility. To this end, 
we make use of a path which is symmetric with the 
pumping step, and consists in two stimulated and one 
spontaneous emissions. To be more explicit, we apply 
two left circularly polarized lasers (see FigHJ) slightly 
detuned from the transitions (60/ < — > 5d,j = |) and 
(5d ; j = | < — > 5p,j = |). Due to these laser fields, the 
atom is likely to fall towards the ground state and emit 
two stimulated and one spontaneous photons. 

Using the selection rules, one can infer that, if a circu- 
larly right-polarized spontaneous photon is emitted, the 
only states to be coupled to the ground level are 
and I72) to and respectively (see Fig©. This 
means that the emission of a right polarized spontaneous 
photon brings the "correct" part of the state vector back 
into Tli — Span^vi) , On the contrary, the other 

cases - "left polarized", "linearly polarized spontaneous 
photon", or "no photon at all"- do not lead to the right 
projection process. 

The "left-polarized photon" and "no photon emitted" 
cases are quite unlikely: indeed the probability that they 
occur is proportional to the square of the error ampli- 
tude, that is to the square of the Zcno interval T, which is 
very short. The "linearly polarized photon" case is quite 
annoying because it mixes the two paths I71) — > 
and j 72) — ► |^2)- This parasitic process and its relative 
probability must be suppressed, with respect to the pro- 
cess followed by the "right -polarized" photon emission. 
This can be done by launching the 78 Rb atom, previously 
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FIG. 4: Coding step through the non-holonomic control technique. The two Hamiltonians H a and H b are alternately applied 
to the system during pulses of timings {U(ns)} ={3.9763, 6.4748, 4.2274, 3.6259, 2.8717, 3.6281, 7.2263, 6.4260, 4.8070, 5.0394, 
6.5242, 4.8890, 4.2400, 7.3834, 4.8653, 5.4799, 4.5341, 4.3099, 6.2959, 3.7346, 6.5293, 6.8586, 6.0749, 5.1213, 4.6806, 3.4985, 
3.9909, 4.6701, 4.5168, 6.4702, 4.7787, 5.3476, 3.4567, 3.8009}. The frequencies of the laser fields involved in the encoding step 
are represented on the spectrum of the Rubidium atom. The fine structure of the Rydberg level 60/ is not represented. 



cooled, into a Fabry-Perot cavity, in an atomic fountain 
manner (fine tuning of the lasers driving the 60/ — 5d 
and 5d — 5p transition will be necessary to avoid reflec- 



tion of the external laser radiation from the cavity). The 
decay rate for the 3-photon transition — ► \i>i) is 



J 



7i"< 



2tt 



7i(Ai + A 2 ) 



r 



2T:hck s 



It 



(y s )., 



where k s is the wave vector of the spontaneously emitted 
photon, 7? r is the left-polarized photon polarization unit 
vector, q (le is the density of states (normalized to 

the cavity volume) for the cavity field at lt s , and the 



bar denotes averaging over the directions of k 8 . The 
transition dipole moments are denoted by d a i,: during 
the projective process the states coupled to I71) and I72) 
are, respectively, 



J 



|Ai) 



5d,j = -,m,j = -- ) ,|A 2 ) 



lMi> 



5d,j = 



5p, j = -zi m 3 = 



|M2> 

I — 



The enhancement (by the presence of cavity) of the the z-axis ensures that 
density of states for the modes propagating paraxially to 

"i ' -^72^2 



d\ 3 ^ k E 2 



ft(Ai+A 2 ) 
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FIG. 5: Decoding step by the non-holonomic control technique. We reverse the magnetic field and the detunings of electric 
fields, as represented on the spectrum of the Rubidium atom, and apply the same control sequence as for coding (same timings) 
backwards. The fine structure of the level 60/ is not represented. 




FIG. 6: Projection path. The lasers involved are marked by solid arrows, the spontaneous photon is represented by a dashed 
arrow. The different polarizations are specified. The fine structure of the level 60/ is not represented. 
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where 7 is the decay rate of \hp,j = §,Wj = +5) into 
1 5s, j = |,mj = +5), so that the undesired process fol- 
lowed by the 7r-photon emission is relatively less impor- 
tant than it were in free space. For the density matrix 
elements p a b the following system of equations can be 
written (i = 0, 1): 

P11I2 — ~ oV'^'Tl 1 '! "^-^72^2)^71721 



V^Ti ^1 -^72 ^2 P71 



72 • 



To avoid dephasing which would corrupt the information, 
the coherence matrix element p 7l72 must be transferred 
with the maximum efficiency into p vlU2 '. the efficiency 



2VT 



71^1 ^72 v 2 



"71"! "I" 1^72 V2 



is thus crucial. According to the Wigner-Eckart theorem, 



C. 



5/2 -3/2 



3/2 -1/2 



C 



3/2 1/2 



3/2 -1/2 1 -1^3/2 1/2 1 -1^1/2 -1/2 1 1 



c, 



5/2 -1/2 



a 



3/2 1/2 



3/2 1/2 1 -l w 3/2 3/2 1 



3/2 3/2 
1/2 1/2 1 1 



where in the right-hand-side the ratio of the products 
of the Clebsch-Gordan coefficients corresponding to the 
transitions stands. These coefficients, which can be found 
in [26|, lead to 77 = 12V2/17 w 0.99827. In other words, 
the probability of error during the Zeno projection stage 
due to the small difference of the Clebsch-Gordan co- 
efficient products for the two paths is equal to or less 
than 1 — r\ w 0.00173 (the equality is reached if the ini- 
tial state is (|0) ± |1))/a/2). Note that the states 60/, 
5d, and 5p have finite lifetimes Tk (see FigEJ - Thus the 
transition rates r 7i „ ; must be much larger than l/Vgo/, 



/r 5d , and 



£'1 



ft(Ai+A 2 ) 



/t5 P , in order 



to diminish errors caused by the decay of these unstable 
states. 

To complete the projection step, one has to transfer 
the atom in its coherent superposition back to the 60/ 
state: this is achieved by the same pumping sequence as 
in the first step. The mismatch of the Clebsch-Gordan 
coefficient products will cause again the error probability 
1 — 77. The information is then restored with very high 
probability and the system is ready to undergo a new 
protection cycle. 

From the beginning of this section we have neglected 
the fine structure splitting of the level 60/, which is 



approximately 2.10 



1 and corresponds to a period 



Tf ~ 1.5/zs. To conclude this section, let us now take it 
into account and see its effect on each step of our scheme. 

Obviously the pumping and projection steps will not 
be affected by the fine structure, since the information- 
carrying vectors {I71) , I72)} belong to the same multiplet 
(J = 5/2). 



The coding and decoding steps are neither modified 
by the existence of the fine structure. Indeed, since the 
typical period of the fine structure Hamiltonian t/ ~ 
1.5/xs is more than 10 times longer than the total duration 
of the coding or decoding steps, it is legitimate to neglect 
its effect. 

The influence of the fine structure on the free evo- 
lution period during which errors are likely to occur is 
more complicated to study in the general case. Yet, two 
simple limiting regimes can be considered. If the spec- 
trum of the coupling functions / m (i)'s is very narrow 
(i.e. if the variation timescale of the / m (t)'s is much 
longer than Tf), one can show that our scheme applies 
directly as though there were no fine structure, provided 
the error Hamiltonians ^£7 m | are replaced by 

where E^' is obtained from E m by simply setting to zero 
the rectangular submatrices which couple the two multi- 
plets (J = 5/2,7/2). The second limiting regime corre- 
sponds to a very broad spectrum for the / m (i)'s (varia- 
tion timescale much shorter than 77): in that case, one 
can show that our scheme applies provided one chooses 
a Zeno interval T multiple of r/. 

In all this section, we implicitly supposed that the Rub- 
dium atom was alone ; but in actual experiments, one 
usually works with an ensemble of atoms : this generates 
serious experimental drawbacks which we deal with now. 
Rydberg atoms are sensitive to Doppler effect : never- 
theless, in the case of cold atoms, this is negligible. But 
the most dramatic effect is due to interactions between 
atoms such as dipolar forces • In a standard Magneto 
Optical Trap containing about 1000 atoms in Rydberg 
states (n ~ 60), the typical energy of these interactions 
is 1 MHz, corresponding to a dephasing time of 1ms [28| . 
As different atoms see different environments, and are 
thus subject to different interactions, it will be impos- 
sible to properly code and thus protect the information 
stored in the different atoms of the ensemble. Beyond 
these problems, we nevertheless want to emphasize the 
demonstrative value of our example : the system consid- 
ered here (Rubidium isotope in a Rydberg state) , though 
not completely satisfactory from an experimental point 
of view, is indeed quite practical for a straightforward 
demonstration of our scheme, since information carrying 
subsystem and ancilla are clearly identified, and every 
step is "simply" achieved. Other systems must be found, 
which will be addressed in future publications ; however, 
the application considered here has already suggested the 
physical relevance and applicability of our method. 



V. CONCLUSIONS 

In this paper, an original scheme has been presented 
which allows to protect the quantum coherence stored in 
a information system X against the action of a set of M 
given error- inducing Hamiltonians E^. 

The information initially stored in the Hilbert space 
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Hi of the information system is transferred into a sub- 
space C of the Hilbert space H = Hi ® Ha of the com- 
pound system T <g> A formed through adding an auxil- 
iary system A called ancilla to the main system. A mul- 
tidimensional generalization of the QZE has been pre- 
sented which makes it possible to protect such a subspace 
against the action of the E^s, provided the dimension A 
of the ancilla meets the Hamming bound A > M + 1. 
The information is thus encoded in another subspace C, 
called the "code space", through the application of the 
coding matrix C: in this appropriate subspace, the error- 
inducing Hamiltonians Ek act orthogonally. After a short 
time, the information thus contains a small orthogonal 
erroneous component due to the action of the Ek ; it 
is then decoded by application of C~ x and restored by 
an appropriate physical measurement which projects the 
state vector onto C with very high probability. The rep- 
etition of this sequence as long as needed protects the 
information stored in the system. 

A physical achievement of the coding and decod- 
ing steps have been proposed which employs the non- 
holonomic control technique. The different algorithmic 
tools needed to implement our scheme have been pre- 
sented. 

Finally, an application has been proposed which makes 
use of the Rubidium atom. One qubit of information is 
encoded in the spin states of the atom whereas the or- 
bital part plays the role of the ancilla. A realistic physi- 
cal setting has been proposed: in particular, a projection 
process based on the spontaneous emission has been sug- 
gested. 
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APPENDIX: EXPLICIT DERIVATION OF THE 
CODE SUBSPACE 

In this appendix we deal with a particular physical 
situation in which the code subspace C can be explicitly 
derived. We consider an atom with zero nuclear spin on 
the level characterized by the orbital angular momentum 
L. The electronic spin of the atom is S = 1/2. The 
natural basis wave functions are \L, Ml; S, Ms). A qubit 
of information is encoded on the two states 



|7i) = \JLSM. h ) 

= E cIm^sms^MliS^Ms), 



1,2, 



where C^m^sms is the Clebsch-Gordan coefficient, 
the scheme we proposed for a Rubidium atom, L = 
J = 5/2, M, ll = -3/2, M. h = -1/2, that is 



In 

3, 



l7i} = 
I72) = 



60/, j = -, mj =-- 



finf . 5 1 
60/, j = -,mj = -- 



We want to protect this information against the action 
of 6 independent error-inducing Hamiltonians Ek, 3 mag- 
netic and 3 electric interaction Hamiltonians. We shall 
see that the code space C = Span [I71) , I72)] can be sim- 
ply built from physical considerations on the action of 
the Hamiltonians Ek- 

Let us first consider magnetic errors. The interaction 
Hamiltonian of the atom with the magnetic field 7? di- 
rected along the fc-th axis (k — x, y, z) is 



/j B B k (Lk + 2S k ), 



u B being the Bohr magneton. Remembering that L a 



L++L. 



2 ' Ly — 2i 

creasing (lowering) the z-projection of the orbital angular 
momentum, and the similar relations for the spin opera- 
tors, one can conclude that a pair of the states with def- 
inite ^-projections of orbital and spin angular momenta 
is a good basis for the code subspace if the difference of 
the of their quantum number Ml = L z is greater than 
or equal to 2. To use this option, one needs to consider 
the error caused by a magnetic field oriented along z. 
The states with definite L z , S z are the eigenstates of the 
Hamiltonian E%. A general superposition of two such 
states must not be rotated in the Hilbert space under 
the action of E% . This means that the eigenvalues must 
be equal to each other. Thus the states 



L+ - L - where L + (£_) is the operator in- 



7f> = \L,M L ;S,M S = +1/2), 



t!) = \L, M L + 2; S, M S = -1/2). 



with M < L — 2. constitute a good code basis for pro- 
tecting one qubit against the action of the E%. 

We shall now consider errors caused by quasistatic elec- 
tric fields. The static Stark shift of a level with zero fine 
splitting, i.e., a highly excited Rydberg level, like 60/ in 
Rubidium, caused by the electric field £ oriented along 
z, is given by const — b£ 2 M^. The value b characterizes 
the polarizability of the atom in the given state. Omit- 
ting the irrelevant constant, we may represent the Hamil- 
tonian of the atom-field interaction (with respect to the 
particular manifold of sublevels of the given atomic state) 
by the operator 



El 



-b£%Ll, k = x,y,z. 



(A.l) 



Ml, Ms 



Note, that since the fine splitting is zero, the spin vari- 
ables are unaffected by the Stark effect. Rewriting the 
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operators E e x in terms of L z , 



errors, too. The basis vectors of this subspace are 



K 

El 



-bSi 



-bet 



4 L+ 



1 



l L ~ + 2 L{L 



1- 



-L(L+1)--L 



7i' E ) = |i ) Af i = -l;S,Af s = +l/2) ! 
= \L, M L = +1; S, M s = -1/2). 



one can conclude that the basis of the coding space can 
be a pair of states of opposite S z and opposite L z (so 
that \L Z \ is the same for both of these states). Indeed, 
these states are not mixed by the Hamiltonian QA.10 ). 
which does not cause spin flips. The error vector is always 
orthogonal to any their superposition, as can be easily 
seen. Among various code subspaces protecting against 
electric errors there is one that protects against magnetic 



It may happen that for singlet electronic states of 
atoms with non-zero nuclear spin, whose nuclear mag- 
netic moment [i n ud is incommensurate with [Ib , one can- 
not apply this explicit derivation of the code space for the 
correction of both the "electric" and "magnetic" errors. 
One has then to look for a more complex coding trans- 
formation through the algorithm we presented in Sec. III. 
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